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(57) Abstract 

A magnetic resonance (MR) imaging ap- 
paratus and technique exploits spatial infor- 
mation inherent in a surface coil array to in- 
crease MR image acquisition speed, resolution 

and/or field of view. Partial signals are ac- ! 
quired simultaneously in the component coils J 
of the array and formed into two or more sig- 
nals corresponding to orthogonal spatial rep- 
resentations. In a Fourier embodiment, lines 
of the Jk-space matrix required for image pro- 
duction are formed using a set of separate, 
preferably linear combinations of the compo- 
nent coil signals to substitute for spatial mod- 
ulations normally produced by phase encoding 
gradients. The signal combining may proceed 
in a parallel or flow-through fashion, or as 
post-processing, which in either case reduces 
the need for time-consuming gradient switch- 
ing and expensive fast magnet arrangements. In 
the post-processing approach, stored signals are 
combined after the fact to yield the full data ma- 
trix. In the flow-through approach, a plug-in 
unit consisting of a coil array with an on board 
processor outputs two or more sets of combined 
spatial signals for each spin conditioning cy- 
cle, each directly corresponding to a distinct 
line in k-space. This partially parallel imag- 

. ing strategy, dubbed Simultaneous Acquisition , . 

of Spatial Harmonics (SMASH), is readily integrated with many existing fast imaging sequences, yielding multiplicative time savings with- 
out a significant sacrifice in spaual resolution or signal-to-noise ratio. An experimental system achieved two-fold improvement in image 
acquisition time with a prototype three-coil array, and larger factors are achievable with other coil arrangements. 
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Simultaneous Acquisition of Spatial Harmonics (SMASH): 
Ultra-Fast Imaging with Radiofrequency Coil Arrays 

Background 

5 This invention is generally related to magnetic resonance imaging (MRI) using 

nuclear magnetic resonance (NMR) phenomena. It is particularly directed to a method 
and corresponding apparatus for more efficiently capturing and providing MRI data 
suitable for use in multi-dimensional imaging processes. 

MRI is by now a widely accepted, medically important and commercially viable 
1 0 technique for obtaining digitized video images representative of internal body tissue and 
structures. There are many commercially available approaches and there have been 
numerous publications describing these and other approaches to MRI. Many of these 
use multi-dimensional Fourier transformation techniques which are now well-known to 
those skilled in this art. 

15 In general, MRI devices establish a constant homogeneous magnetic field, apply 

a specific additional bias field gradient in a known plane or region under consideration 
to orient nuclear spins, and apply a radiofrequency pulse or a sequence of pulses to 
perturb the nuclei. These nuclei in the known bias field gradient emit an RF signal in a 
specific band determined by the magnetic field distribution, and these RF emissions are 

20 defected by receiving coils and stored as a line of information in a data matrix known as • 
the k-space matrix. 'The full matrix is built up'by successive cycles of conditioning the 
*spins, perturbing them,- and collecting RF emissions. An image is then generated from 
this matrix by Fourier transformation, which converts the frequency information present . 
ih the RF, oscillations to spatial information representing the distribution of nuclear spins;; 

25 in tissue or other imaged material. . ^ ' ' ; ' 

* * Magnetic resonance imaging has proven to be a valuable clinical diagnostic tool 
in a wide* range of organ systems and pathophysiologic processes. Both anatomic and 
functional information can be gleaned from the MR data/and new applications continue 
to develop with each improvement in basic imaging technique and technology. As 
30 technologic advances have ^ improved achievable spatial resolution, for example, 
increasingly finer anatomic details have been able to be imaged and evaluated using MR. 
At the same time, fast imaging sequences have reduced imaging times to such an extent 
that many moving structures can now be visualized without significant motion artifacts. 
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Often, however, there is a tradeoff between spatial resolution and imaging time, 
4 since higher resolution images require a longer acquisition time. This'-balance between 
^ spatial: and temporal resolution is particularly important in cardiac MR, where fine 
- details of coronary artery" anatomy/ for example,* miist be discerned on the Surface of a 
5* rapidly beating heart, -A high-resolution image acquired over a large fraction of the 
: cardiac cycle will be blurred and distorted by bulk cardiac motion, whereas a very fast 
image may not have the resolution necessary td trace the course and patency of coronary 
arteries. Some; of the fastest imaging sequences- currently implemented, such as echo 
planar imaging (EPI), approach the goal of yielding images of reasonable resolution in a 
10 suitably short fraction of the cardiac cycle. Other approaches have also been tried to 
eliminate-thte effects of cardiac motion; including £-space segmentation, in which image 
acquisition is divided up over several cardiac cycles with ECG gating to ensure that the 
heart is in the. same phase of systole or diastole : during acquisition of each segment. 
: Cine images, of multiple cardiac phases may be pieced together with this technique, with 
.15; . partial acquisitions for different phases occurring in each cardiac cycle. One problem 
with this class of techniques is that respiratory motion can change the position of the 
, , .heart over the course of several cardiac cycles. Partial acquisitions will then be 
misregistered, and. artifacts will result In A an ; attempt to . eliminate or adjust for 
• respiratory motion, 1( breath holds, respiratory gating, .and navigator . echo gating 
, ,20 i; . techniques have all been tried, and each, of, these techniques, has had spme significant 
. . ( successes, Nevertheless, .an imaging strategy which allowed high-resolution images to 
be acquired comfortably withi^ one or two phases of the cardiac-cycle would circumvent 
^ f i: many^ of . the difficulties and residuah . artifacts associated with these compensation 
: , ^ techniques. , ...... . .. , t> . , 

25 Th e speed with 'which magnetic resonance (MR) images may be : acquired has 

I - v.: ^Jready- increased dramaticaljy ; oyer the^past decade. The improvements in speed may be 
^ traced to a- combination of advances in- the technologies of magnet construction and 

: i .-actuation, and innovation? in imaging strategy. Strong, fast-switching magnetic field 

gradients and fast electronics have allowed the intervals between data collections to be 
- ,30. , reduced significantly. Meanwhile,, fast gradient-echo and spin-echo sequences, have 
reduced image , acquisition time by allowing ^greater portions , of A:-space to be sampled 
after each spin excitation. Echo planar imaging (EPI), fast low-angle shot (FLASH), 
turbo spin echo (TSE), and spiral imaging techniques alJ allow very short intervals 
between acquisition of successive data points. The DUFIS, OUFIS, RUFIS, and 
35 BURST family of sequences further reduce image acquisition time by eliminating time 
delays incurred during gradient switching and echo formation. Details of the above- 
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menlipned eight techniques may be found in the following papers: P. Mansfield, Multi- 
planar image formation using NMR spin echoes . , J: Phys. C. ,10, L55-58 (1977); A. 
Haase, J. Frahnr, D r Mattaei, ,W. Hanjcke r K. D. Merboldt FLASH imaging: rapid NMR 
imaging using low flip-angle, pulses . J. Magn. Resort. 67, 256-266 (1986); J. L. Listerud, 
, 5 . S. Einstein, E. .Outwater, H., Y. Kresseb .First principles of fast<spin echo. t Magn. Resort. 
Q. 8, ,1 99-244. ,(1992);. G. e Meyer, ; B, ? Hu,,<D: : Nishimura, . A. Macovski, Fast spiral 
coronary artery , imaging.^ .. Magn. r Resort. Med: 28, 202-213 (19923; I ; J> Lowe, R. E. 
Wysong v DANTE: ultrafast imaging sequence (DUFISV J. Mazn: Respn. Ser. B 101, 
r 106-L09 (1993); L. Zha, .1. J. Lowe, Optimized ultra-fast imaging sequence (OUFIS). 
10 Magn. Resort. Med. 33, 377-395 (1995); D. ; P. Madio, I. J. Lowe, Ultra-fast imaging 
, using low flip angles and FIDs.. Magn.. Reson. Med 34, 525-529 (1995); and J. Hennig, 
M. Hodapp, Burst imaging . MAGMA 1. 39-48 (1993V 

■ " * Increasing the speed of MR imaging further is a challenging proposition, since 

the aforementioned fast imaging techniques' have already achieved an impressive 
15 efficiency. All these techniques allow very short intervals between acquisition of 
successive data points, and hence do' not waste* much' time in accumulating the data for 
the ^-space matrix required to generate an image. Iri flow-encoded EPI images, for 
example, the entire complex £-space matrix is filled in a single spin' excitation (which is 
' followed by multiple spin conditioning cycles involving the application" of multiple 
20 / stepped field gradients), and' the resulting image matrix is likewise "full," with useful 
information stored in 1 both the real and the imaginary channels.' One common feature of 
hearly ail the fast imaging techniques curreritly iri use, however, is that they all acquire 
-data in a sequential fashion^ - Whether the required data set, i.e.,' the k- space data matrix, 
is filled in a rectangular raster pattern, a spiral pattern, a rapid series of line scans, or 
25 some other novel .order, it is acquired one point and one line; at a time. r y : ; 

J 11 ; That is, the prior aft in fast MR imaging has' coftcentrateel on increasing the speed 
of sequential acquisition by reducing' the internals between scanned lines.' Modifications 
to pulse sequences or to magnetic field gradients have yielded a gradual improvement in 
imaging y - speed by allowing faster sequential 'scanning of fc-spacfe,; but these 
30 improvements face limits due to' the intervals necessary to create, switch or measure the 
'magnetic fields or signals involved in data acquisition. It would therefore appear 
■ ; difficult to devise & sequential technique with : significantly better efficiency than the 
-current fast imagingtechniques. 

Several fast imaging schemes have been proposed to . date using simultaneous 
35 * data acquisition in multiple RF coils, as described in: D. Kwiat, S. Einav, G. Navon, A 
decouplecl coil detector array for fast image acquisition in magnetic resonance imaging . 



WO 98/21600 PCT7US97/20068 . 

& 

-4- 

MedPhys, 18:251-265 (1991); D, Kwiaty S. Eihav, Preliminary experimental' evaluation 
of an inverse source imaging procedure- using a decoupled -coil detector array in 
"magnetic resonance imapinp, ' - Med Rv>d Phy, - "j t 257-263 (1995); J; W. Carlson, T. 
Minemura, Imaging time reduction through multiple receiver coil data acquisition and 

5 > ima S e reconstruction. Magn ^7?e^/?-Me^-29r681-688 (I9'93) and^ U.S. Patent 4,857,846; 
arid J. B. Ra, C. Y. Rim, F£st imag inq f; ubirig l subenb6dihg data sets from multip le 

' ^z\ohrMagH%eson htiedM, 142- : 1 45 ( 19^3):-^ ThesS approaches have offered the 
prorniseof significant savings in Tm^ -'t-'-^' : 

Carlson and Minemura describe a twofold acquisition time savings using 
10 two nested body coils. "In their approach, partial" "data sets are" collected 
simultaneously in the two coils, one of homogeneous sensitivity and the other 
with a linear gradient in sensitivity.Missing lines in A>space are generated using a 
series expansion in terms of other phase-encoded lines. This approach using 
body coils appears to require that a significant portion of the data for the partial 
15- k-space matrix be acquired before any of .the missing lines can be filled in by - 
t -.h postprocessing, and thus does not allow for the missing lines to be built up as the 
'•_ .data arrives,, in real time ; The approach uses coil sensitivity information in place 
, of : some portion of the gradient phase encoding steps, but has drawbacks. The 
r coils used by. Carlson and Minemura are body coils, which provide. large volume 
20 coverage but lower overall. sensitivity, than.surface. <x>Us,.-and it wpuld be difficult , \ 
> : to augment their number to improye,time savings..- ;: 

THfe approach of Ra L and Rim involves a simultaneous acquisition technique in 
"'which images of reduced FOV are Acquired in multiple coils of an array and the Nyquist 
' ■ aliasing in those images is*' undone" by -reference ^ to component 'coil sensitivity 
'"^S'^irifomationV- Tlle'unaiias r irig T "procedure 'involves a pixel-by-pixel'' matrix inversion to 
Y regenerate the ; full FbV' from multiple copies of the aliased image data. The 
'~ ; " "subencoding" tebhriique 3 6Y'R* : and Rim * relies ; ori estimates' bf component coil 
sensitivities by effectively probing 4 the' sensitivity at each pixel. This pixel-by-pixel 
'[ approach can lead to local artifacts; for example, the matrix inversion can begin to fail in 
30 regions of loW sensitivity. Further, by its Very nature as a pixel by pixel ' dealiasing 
approach, the Ra & Rim method is computation-intensive and 'ik limited to 
postprocessing, as all image data must be present before the reconstruction can be 
undertaken. - ~ * - ! - - 1 / ; 

In a related area, multiple coil signal collection has been used in MR phased array 
35 systems as reported in R.B. Roemer, W.A. Edelstein, C.E^ Hayes; S.P. Souza, and O.M. 
Mueller, The NMR phased array. Magn. Reson. Med 16, 192-225 (1990); C.E. Hayes 
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and P.B. Roemer. Noise correlations in data simultaneously acquired from multiple 
surface coil arrays . Mazn. Reson. Med 16, .181-191 (1990); C.E. Hayes N. Hattes, and 
P.B. Roemer-j Volume imaging with MR phased arrays . Magn. Reson. Med. 18, 309-319 
.(1991). -The. increased information content of the multiple received signals in such 
5 ^systems has been used to increase the signal-to-noise. ratio (SNR) of MR images. Since 
, their initial description, phased. aijays have; seen increasing use in. clinical MR imaging. 
• For example, .the improvements An SNR, provided by phased arrays .have allowed 
significant advances in imaging, of t the pulmonary vasculature as reported. by T. K. F. 
Foo, J. R. MacFall, C. E. Hayes, H. D. Sostman, and B. E. Slayman, Pulmonary 
10 vasculature: single breath-hold MR' imaging with phased array coils . Radiology 183, 
473-477 (1992). Still, with a few 'notable exceptions, the bulk of phased array 
applications have addressed increased sensitivity, with little effort towards improving 
image acquisition speed or resolution. 

Summary of the Invention , . 

15 " An MRI system according to the present invention uses a multiple-coil data 

-collection system to acquire some portion of the k-space matrix in parallel, rather than 
sequentially in tibie. ; In a preferred embodiment; signals are obtained from multiple RF 
coils each occupying a different position -with respect to the imaged volume, and each 
therefore having different but at -le'ast partially overlapping spatial sensitivities. The 
20 multiple coils are preferably positioned and/or -their outputs sampled in a manner to 
minimize inductive coupling; . ,but -they need riot individually - span the - full region of 
f .. interest nor be fully independent. The. signals collected in this plurality of coils are then 
combined with suitably chosen weights, to produce two on more composite signals, each 
. , of which approximates a sinusoidal sensitivity, modulation,, or . "spatial harmonic." As 
.25. used herein, the. term "spatial harmonic" refers to ; a sinusoidal and/or. cpsinusoidal 
, ^variation in spatial, sensitivity witlva wavelength. that is an integer fraction^of the extent 
of the field of view. Each line of these spatial harmonic composite signals constitutes an 
< additional line of the £-space matrix which would require a distinct gradient step in a 
.,. , conventional MR. acquisition. Thus, the use of such signal combinations eliminates the 
, 30, need for some of the conventional gradient steps. The technique just outlined is referred 
7 to herein as SiMultaneous Acquisition of Spatial Harmonics (SMASH), and it may be 
. used to reduce image acquisition times by a multiplicative factor without a significant 
sacrifice in spatial resolution or signaUto-noise ratio (SNR). The invention also includes 
non-Fourier embodiments, in which the coil signals are tranformed or combined with 
35 weights to yield composite signals which each correspond to a non-Fourier spatial 
distribution, such as a wavelet". 
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The SMASH technique in its prefeiTed embodiment is implemented' in an MR 
imaging apparatus and uses linear combinations, of- the simultaneously acquired signals 
, . • from multiple surface coils to .generate multiple data sets with distinct offsets in £-space. 

The full A-space matrix is reconstituted, by r interleaYing;these ? shifted. data sets, and an 
5 image is therefore generated with only a fraction of the usual number of gradient phase 
encoding steps. Consequently,' the total image acquisition i time may be reduced or the 
total quantity of data generated during a fixed acquisition time may be increased by this 
same fraction- * . ••'* j ..-r - , • .' " 

Conceptually, 'the* invention may bV ; se'eh a ; s partially replacing gradient phase 
.10 encoding by a . spatial encoding,, procedure tied to the detection apparatus. In SMASH, 
some of the spatial modulations that normally distinguish differently phase-encoded 
lines are generated instead by amplitude modulation, resulting from combinations of 
component signals from multiple coils arrayed above or around the imaging volume. 
Computationally, 'by shifting responsibility for spatial encoding from the spin 
15 preparation stage to the stage of -signal detection and combination, applicant is able to 
combine the multiple simultaneously acquired signals, either immediately 6h-the-fly or 
after-the-fact, .to , produce multiple different spatial harmonics.- In this manner, multiple 
lines in it-space are acquired simultaneously, in a parallel rather than a purely sequential 
acquisition scheme. 

•.•-•7-?P; -u v: .'• A representative apparatus irnplementing.thejnyention simultaneously acquires 
partial signals from multiple coils in a surface coil array and combines them into two or 
more differently weighted combinations.that accurately represent several spatial . : - 
harmonics. The weights may be generated by theoretical calculations based on coil 
geometry, or else they may be derived by a calibration protocol carried out on a phantom 
c25 : - image br'bri a pre-scan at the time of in vivo imaging! The calibration protocol uses a 

numerical optimization algorithm to determine coefficients for linear signal combination 
which will best approximate the desired spatial harmonics. £>nce the spatial harmonic 
signal combinations have been formed, the missing lines in k-space are filled in, and the 
conventional MR image is generated by Fourier transformation of the expanded data 
30 matrix. Preferably, the coils in a SMASH coil' array are. surface coils, i.e. coils 

positioned on the surface of the body and designed to capture the MR signal efficiently 
over a restricted region of interest. Such surface coils typically have a spatial extent on 
the order of 5-50 centimeters. 
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Brief Description of The Drawings - 

' ' - These and other features of the 'invention* will' be understood from the description 
herein read in 1 light of the prior art, together with the drawings illustrating representative 
-embodiments of the- invention, wherein - : 

5 FIGURE 1 illustrates a general MRI apparatus having features common to the 

invention and to the' prior art; 

FIGURE 2 illustrates a conventional data acquisition scheme; , 

FIGURE 2 A. iljustrates data acquisition of the present invention 

FIGURE 2B illustrates a receiver coil array of the ihvehtibn and corresponding 
10. image- plane geometry; - • . - ' * • ; 

FIGURE 3 illustrates' another receiver coil array of the invention togetherwith 
coil sensitivities and the construction of spatial sensitivies from component coils; 

FIGURE 4 illustrates images formed from the coil array of FIGURE 2B and 
4 _ complex sensitivities ofthe coils; ....... . .. : . . ( 

15 FIGURE 4A illustrates fitting spatial harmonics with 1 weighted coil sensitivity 

f ■- - . functions; . r- : ■ *: . . . . - r \- • . >•* • •* ■ - 

FIGURES 5 and 5A show processing and representative images with the signal 
1 Trom a set of non-optimized coils and the fitting procedure of FIGURE 4A in a prototype 
- . embodiment of the invention;** - : : " ■"*'■<■- . : ' : • 

20 FIGURES 6A and 6B show MR images of a phantom made by processes of the 

prior artand' the invention, respectively ; and " "/ 

^ r FIGURES 7 A and 7B. show. MR images of a person's head^mad^ by process 

. ( the prior artand the invention, respectively; : ... £ - : , - ^ ( -\- 

> ■•■ . Detailed Description « J / : i ■ 

25 ' ' ' FIGURE 1 illustrates schematically an MRI system 10 which includes 

the usual static magnet assembly, gradient coils, and transmit RF coils 
collectively denoted 12 under control of a processor 14, which typically 
communicates with an operator via a conventional keyboard/control workstation 
1 16. These devices generally employ a system of multiple processors for 

30 carrying out specialized timing and other functions in the MRI system 10 as will 
be appreciated. Accordingly, as depicted in FIGURE 1, an MRI image processor 
1 8 receives digitized data representing RF NMR responses from an object region 
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15 



20 - 



under examination (e^g., a human- body 1) and, typically via multiple Fourier 
transformation processes well-known in the art, calculates a digitized visual 
image (e.g., a two-diniensional array of picture elements or pixels, 1 each of which 
^ may have different gradations of gray values tit color values, of the like) vvhifch is 
' then be converitiohal5y : displdyed, or printed 1 out, : on a display 18a. - ? v 

In such overall operation, the apparatus of the present invention is 
largely conventional. However, in accordance with the present invention the ' 
basic RF : data acquisition is modified, and subsequent signal processing altered, 
* by providing f a plurality r t>f surface coils 20a/20b..:20i for simultaneous signal 
reception, along with corresponding signal processing arid digitizing 1 channels. 
The processor recbmbines the collected values into two or more spatial 
haririonics fromwhich multiple lines of the signal matrix are developed. This 
recombination may tie performed in real time as the data arrives, or after th^ fact 
via postprocessing as is convenient with the apparatus and the calibration 
information at hand. ' ; ; ' 1 - ; - 

Operation of the invention will be better understood from a description of 
gradient phase encoding in conventional processing and a comparison to the 
spatial encoding' produced by amplitude modulation of the multiple coil signals 
' in the present invention, and a consideration of the inforrnation available 1 in 
spatial harmonics. : * h ^ : 2 5 v 

In the g eneral case, the ma gnetic re sonance signal for a plane with spin 



density |P( x >.y) |and coil sensitivity l C ^ x ^1 may be written as follows: 



25 



30 



|s(k x ,k y ) = ffdxdyC(x,y)p(^:y)e X p{-ik^x-ik y y}= Cj^kJ® p(K\) 



[1] 



where 
and 



as ..usual, , with g the gyromagnetic ratio, 
Jthe magnitude of the x and y gradients, and es spent the 

in the x and y gradients, respectively. Here, the spin excitation function as well 
as the effects of relaxation have been incorporated into a pulse-sequence-specific 
sensitivity function C. The symbol EE] indicates a convolution. For regions of 
the sample in w hich the c oil sensitivi ty is roug hly homogeneous we may write 
arid 



C(x,y)= 1 



S(k x ,k v ) 



is equal to 



p(k*,M 



the spin density function. Double Fourier transformation 
TO « 
reconstructs the usual spin-density image 



, the spatial Fourier transform of 

Ik I * 

with respect to L*J and 



P(x,y) 
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, ( In a prior art phase-warp imaging method, a k-space,. matrix of frequency- 
encoded and phase-encoded information, as represented generally in FIGURE 2 
\ is produced for the image plane of interest, with successive lines represented by 
broad, horizontal line corresponding to a different phase-encoding gradient value. 
5 Each phase encoding step corresponds to .a .distinct spatial, modulation and these - 
spatial modulations are represented schematically by the curves at the left of the 
FIGURE. • < r . . :>j 

■ > N ■ Generally prior, art magnetic resonance receiver :coils, especially, surface 
coils,, do not have uniform sensitivity. Signals from different regions of the 
10 imaged volume produce different currents .in an RF coil, with the spatial , 
variation in sensitivity being simply related, to the inhomogeneity .of RF field 
produced by the coil over the sample volume. For a standard, circular surface 
, .cqil,. there, is a. sensitivity "s\yeet spot" centered at roughly one diameter below 
the coil, with a monptonic falloff of sensitivity together with increasing- phase 
15 differences in all directions. Traditional imaging protocols often position the 
receiving coil with the target tissue at its region of maximum sensitivity. 

Applicant realized, however, that when collecting a signal with some 
combination of coils having an aggregate sinusoidal and cosinusoidal spatial 
sensitivity profile, the MR signal froni these coils will have, an information 
20 content somewhat different from that of the usual coil signal,, and .that by., 
separating out one or more collected signals corresponding to pure spatial 
harmonics,' these may be used l io fill a* larger portion of the "data -.space than is 
conventionally done. 1 * A complex combination of cosinusoidal arid sinusoidal 
sensitivity profiles yields 



.25 



G( x„y) = exp {iK^ y } = cos Ky y ; +i sin K y y 



[2] 



where 



is the spatial frequency of the inhomogeneous coil sensitivity . -For the 



it MR signal, we then have 



S(k x ,k y ) - {Jdxdyp(x,y)exp{-ik x x- i^. - ^ K y )y} = p(k x ,k y - K y ) 



[3] 



In other words, the combined MR signal from the inhomogeneous coils is . , 



■ * , . . . ^ — K 
30 t shifted in /r-space by an amount L d.. This £-space shift ;has precisely ,-the same 



form as the phase-encoding shift produced by evolution in a 3; gradient of 



magnitude 1^ y >' ~ y l -Thus, when this shifted signal is acquired 
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simultaneously with the usual signal from a horriogeneous coil, one gains an 
extra line in &-space for each application of the x and y gradients. Applicant 
exploits this realization to 1 convert the signals collected in si set -of surface coils 
- into spatial harmonic signals, and fill multiple lines of k-space from each signal 
5 collection. If one such harmonic is generated in addition to the homogeneous 
"- signal (zero th harmonic) then one' may then apply^alf as many 1 ^ gradient steps 
to build up the full A-spa£fe matrix necessary for image acquisition. When higher 



harmonics of larger are added, the number of necessary y gradient steps is 
' further reduced. For a multiple coil setup representing Ki spatial harmonics in 
1 0 addition to the homogeneous signal, the acquisition time is reduced by an integer 
factor of M+L If both positive and negative exponentials are used, ■ 

[4] 



C±(X y) = exp {±iKyy } = cos Ky y £ i sin Ky y 



an extra line appears at an intervalL_ll above and below each gradient-encoded 
line. The £-space geometry of this acquisition strategy is illustrated in FIGURE 
15 2A. 

Thus, if one were able to collect a signal with a coil shaped so that the 
response /it picked up had the form of such a spatial harmonic, it would directly 

* ; yield' a k-space entry corresponding to oner of the usual spin-conditioned 

* .:, measurements; 'an array of such coils might simultaneously acquire several lines 

20 of data. ; , , ,. • . ....... ... . • ' • ,' ; j , ... . 

Applicant further realized that the sensitivity of each individual coil in an 
■ airray need not be strictly sinusoidal, so long as the 'net sensitivity of the array 
permits construction of *harmoriics having the j desired ^ sinusoidal 4 shapei This 
"realization significantlyfrelaxes the constraints "on coil design arid 'disposition, 
25 and in conjunction with a calibration or normalization process, allows the signals 
._. ; ,:from multiple coils with a wide range of shapes and geometries to be combined 
. in several ways to yield net signals with multiple sinusoidal sensitivity profiles. 

FIGURE 2B illustrates geometry of the coil array and the image plane. A " 
set of three adjacent coils are used with separate coil outputs in a surface array 
30 extending generally over the region and planes of interest. The patient is 
positioned and the spin preparation fields are applied to condition a plane P 
which, by way of example may intersect .the patient's heart to image the blood 
vessels thereof, or which may intercept a region of the abdomen to image its 
contents. The coil array 20 is located above or below the region to be imaged, so 
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that each coil 20a, 20b, 20c has at least some sensitivity to. RF signals emanating 
from region P. t . In t the illustrated embodiment .the.- coils-. are each somewhat 
overlapped in the y-direction to minimize inductive coupling. As shown, they are 
overlapped ; at neighboring edges.- Each coil, viewed individually has a sensitivity 
5 function .which as ^highest directly, above or below the -center of the coil, and 
which falls off f with distance from the coil center. : In accordance with the present 
invention, the signals from these, coils, are combined to produce; several .separate 
"virtual". or synthetic signals which each correspond to a pure spatial harmonic. 

< FIGURE 3 demonstrates this situation schematically fori a set of eight 
10 rectangular coils 20a, 20b;... 20h laid end-to-end, with a slight overlap. As shown, 
in line (A) of the Figure^each coil 20a, -20b... haSia sensitivity curve a v b,... which 
rises to a broad peak directly under the coil and drops off substantially beyond 
the coil perimeter. The sum of the coil sensitivities forms a relatively constant 
sensitivity, over the full length of the array, corresponding to the zero 1 * 1 spatial 
15 harmonic. The remaining lines (B)-(E) of FIGURE 3 illustrate recombining of 
different ones of these individual offset but otherwise identical coil sensitivity 
functions into a new synthetic sinusoidal spatial sensitivity. Different 
weightings of the .individual component coil sensitivities lead to net sensitivity 
profiles approximating several of the spatial harmonics of Ky : In the FIGURE, 
20; coil sensitivities (modeled schematically for this Figure as Gaussian in. shape) are 
combined . to produce harmonics at .various fractions of the fundamental spatial 



A. *^7t / IC 

wavelength I y d, with L^jJ being on the order of the total coil array extent 1 

iny. Weighted individual coil sensitivity profiles, are depicted aS;thin solid lines 
^ beneath ; each, component coil. Dashed lines represent- the .sinusoidal or 
25 . * cosinusoidal- weighting functipns.v Combined sensitivity prpfi indicted by 

:thick solid lines. These combined profiles, closely approximate ideal spatial 
harmonics across the array. - . : . • - . v\ 

For a weighted sum of component coil signals, the net sensitivity profile 
is a linear combination of the intrinsic sensitivity profiles of the 



X 



G ,01 -(x;y) 



30 component coils: In this case, for an ^/-component coil array, we have 

• [5] 



C' 0 \x,yl=^n j C j {x,y) 



' where nj is the* weighting coefficient of the j*h coil sensitivity function Cj(x,y) 
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In accordance with the aspect of the invention illustrated in FIGURE' 3, 
by simply weighting the signals in the coil array with^ appropriate weights and 
combining them, one obtains a signal whose amplitude is modulated by the 
spatial harmonic for that combination-of-theeoil sensitivities. That is, the N (in 
5 this case, eight) component coil-' signals' are acquired simultaneously and- then 
multiplied by various weights-and recombined in a total of up to W independent 
combinations each representing -a spatial harmonic corresponding- to a distinct 
gradient phase encoding step . or. distinct offset in k-space. f .This additional 
processing does not add in any way to,, image, acquisition time, and it may be 
1 0 performed quickly, on stored signal data as part of a post-processing algorithm or 
else may be implemented in a simultaneous flow-through fashion as described 
below > .. 



20 



In actual practice, the sensitivity "profiles of RF' surface- coils'^are not 
simple Gaussian profiles bufmore complicated functions which are, in general, 
15 complex in the mathematical sense of having ' both real and imaginary 
components? The coil sensitivity functions must describe both the magnitudes 
and the phases of the signals produced by precessihg spins at various distances 
from the- coil center* and these magnitudes and' phases vary according to the 
reciprocity relation - ; : ..• ... 



^ r >^-^{ B !,v(r) : »n(r)} 



[6] 



where -£(ir) is the voltage induced- iff a coil by a given voxel at position r, m(r) is 



the nuclear magnetic moment of the Voxel, ahd l B '»r( r ')| is' the xy vector 
comporierit of the field generated atV by a unit current in'the coil 'as described for 
example in D. I. Hoult and R. E:' Richards', the signal-to-noise ratio of the 
25 nuclear magnetic reso nance experiment J. Magn. Reson 24; 71-85 (1976). 

; .This relationship- has two main : consequences for applicant's SMASH 
reconstruction: first, applicant assigns complex' : weights nj to capture the full 
information content of the coil sensitivity functions, and second, the target 



complex exponential sensitivities 



C±(x,y) = exp{+/7^>} = cos K v y± i sin K y y 



30 are preferably developed directly, rather than by separate generation and 
subsequent combination of sines and cosines. The real and imaginary parts of 
the complex coil sensitivity combina tions may be fitted simultaneously to the 



real a_nd imaginary parts of | C * | , This may be gone,, for example, with a 
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numerical optimization algorithm* using .the., sum of absolute magnitude 
deviations from the target profile or some other similar quantity as a measure of 
goodness of fit. In the prototype described below, f this fitting was performed in a 
straightforward manner and -integrated with data reconstruction algorithms 
5 . written by applicant, in the MAJLAB- ^programming^ language (Mathwdrks, - 
: Isfatick Massachusetts).,^ Applic^Khas.;now confirmed that such L processing 
produces good image quality at ^greatly enhanced. speed: > ; ^ . ^ * 

■ ' " In setting up ah actual MR iriiaging- apparatus, the coil geometry for 
SMASH signal acquisition may be that of an MR phased array, as described in P. 

10 B. Roemer, W. A. Edelsiein, C. E. Hayes, S. P. Souza, and O. M. Mueller, The 
NMR phased array . MagM. Resort. Med. 16, 192-225 (1990); C. E. Hayes and P. 
B. Roemer, Noise correlations in data simultaneously acquired from multiple 
surface coil arrays . Magn. Resort. Med. 16,181-191 (1990); and C. E. Hayes, N. 
Hattes, and P. B. Roemer, Volume imauint; with MR phased* arrays . Magn; 

1 5 Resort. Med. 18, 309-3.1 9 (1991). Indeed, many of the hardware components 
useful for SMASH imaging, are. already present in traditional phased arrays, 
which contain multiple inductively, decoupled coils with some spatial separation, 
and;which include separate receivers for independent collection of data from the 
coils. The technical problem of minimizing inductive coupling of such coils has 

20 been previously addressed for multiple-coil constructions in the context of 

phased-array MR imaging devices, and two 1 basic strategies were developed. One 
strategy is to design a coil array with an appropriately chosen overlap of 
neighboring component coils to minimize. inductive. coupling. The second ^. iy * 
strategy involves the use of low input impedance.preamplifiers on each ■ . ^ . _ ; 

25 component coil channel. Both of these features are advantageously, applied in ... . 
various embodiments of the present invention. . f . : . ^ 

A Prototype and ^Results. : : r -.' : '-L- — *- : -J ; ■-' 

•. »• Applicant has implemented " a- basic embodiment; of the invention as 
outlined above in the following manner. First, a linear RF coil array as shown in 

30 FIGURE 2B was selected having a geometry suitable for spatial harmonic 
. generation. This array was a cardiac imaging array having three eight-inch (200 
mm) rectangular coils adjacent to each other. The appropriateness of the coil 
geometry was first tested in numerical simulations, using analytic integration of 
the Biot-Savart law to calculate the transverse field Bj^,(r) of Eq. [6], and 

35 hence to model' t&e sensitivity profile of each coil. Taking into account the 
number of spatial ^harmonics that could comfortably be jgenerated using the coil 
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airay, and therefore the fraction of the total A:-space matrix to be collected, partial 
image acquisitions were planned, and a regimen of gradient step's and RF pulses 

- were established. Image data 'were their acquired from each boil 1 simultaneously 
in the separate channels of the array, using a conventional fast imaging sequence. 
5 : Since^'only a fraction of the usual signal data was'to be used to generate 1 the final 

' images - only a fraction of the usual imagihg'ume'was spent oil data collection. 

Next, the stored individual coil image d'ata'were exported to a workstation 
for postprocessing. .Individual complex sensitivity profiles were "determined for 
each of the component coils in the array, usihg phantom image, data. instead of 
10 numerical simulations. This process is illustrated In FIGURE 4. For the 
phantom image technique a round water-filled dummy was suspended at the 
image plane used in the image acquisition above to provide a phantom image 
slice known to be of uniform spin density. Line (A) of FIGURE 4 shows the 
images ajb,c generated from the separate signal received in each of the respective 
1 5 coils 20a, 2*0b, 20c of FIGURE 2B, processed separately. Line (B) of FIGURE 4 
' shows the corresponding real and imaginary components of the coil sensitivity 
'function determined from the detected signal, ^takeri along the diameter A of the 
FIGURE. ■•■ ■.:■>:■■> .-. «- ' > . . 



20 



25 



To determine suitable weights for forming the desired spatial harmonic 
sensitivities from these three coils, an iterative fitting procedure was employed to 
represent each spatial harmonic as a linear combination of the three sensitivities. 
Thjs^ was done as shown in FIGURE' 4A. Firsf, the real arid imaginary 
components of the sensj ; tiyjty"^ctipn' were normalized,' yielding the component 
. -Y5 lues . , a ' b 'f P^" ed M pane's A and B7 respectively,' of FIGURE 4 A. Then, for 
„ $f 2'0itfmm- field of view, optiipal weights for linear combination of component 
coil signals were "determined by iterative fitting of the* sensitivity data to the 
target spatial harmonic sensitivity profiles, with the complex weights serving as 
fitting parameters and the sum of absolute' magnitude deviations from the target 
profile serving as a measure of goodness of fit. Panel Cpf FIGURE 4 A shows 
the resulting best fit to the zero* harmonic, corresponding to uniform sensitivity. 
Panel D .illustrates the best fit to the first spatial harmonic.' Each was realized as 
a simple linear combination, with complex coefficients, of the three coil 
sensitivities shown in panels A and B. 

In accordance with this aspect of the invention, the stored coil signals 
35 received during each spin conditioning cycle of image acquisition were 
recombincd with the weights as determined 'above into' two distinct spatial 



30 
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harmonics, corresponding to two distinct £-space 7 shifted isignals. The data sets 
were then. interleaved, in the ordering sequence. of FIGURE 2A , to yield the full 
fc-space matrix: This matrix; was double Fourier. transformed in the conventional 
manner to give the reconstructed SMASH image. In the. prototype, this data was 
5 > reconstructed after thefact, i.e., post processed to. both determine a suitable set of 
weights for forming the. composite v. spatial harmonic sensitivities, and then to 
recombine the already measured signals with those weights into corresponding 
spatial harmonic measurement signals. In various embodiments of the device, 
however, the optimal weightings may be determined in " advance withprior 

10 knowledge of the coil sensitivities, and the conversion of coil signals may be 
performed on-the-fly, that is, in real time. For example, a simple controller 
attached to the coils may determine the appropriate weights, from a lookup table 
or by some other linear or non-linear interpolative mechanism, and then combine 
the signals to produce a separate composite signal corresponding to each spatial 

1 5 harmonic. A neural network may be used to "learn" the appropriate weighting. 

For simple ,coil configurations, .the sensitivity may be modeled with 
. accuracy, rather than be, empirically derived from normalizing the coil responses 
to a uniform sample. More generally, the process illustrated in FIGURES 4 and 
4A may be used to determine suitable coefficients for a given image plane, and 
20 this process may be repeated to compile. a table for each of many, planes so that 
the apparatus need only look up the necessary weights in order to produce the 
desired composite signals for the plane, of interest. When the sensitivity 

1 information is known in advance, then, reconstruction does not require a 

, - ■ -■ * ■ ' * ' v & rr. m « ' ; -'. -\\\ , 

preliminary sensitivity measurement or an iterative fitting procedure, and simply 
- s' . - I > • ' % ; K- *"'.->■ •"■•.* ■■• - r -- '"' -' - ■ • i-: ' V- ' - 

25 involves a set of weighted slims. In either case, once the new composite, signals 

are formed, the additional data must simply be entered at the correct matrix 

position, and the full matrix subjected to ia fast Fourier transform. This 

processing is as straightforward as the postprocessing which is commonly 

performed when conventional images are generated from phased arrays. 

30 Furthermore, once the full k-space matrix is constructed using the SMASH 

technique subsequent processing of the k-space matrix to yield the image is 

identical to that presently employed in substantially all Fourier MR imaging 

systems... 

FIGURES 5 and 5A illustrate the k-space matrix construction for this 
35 prototype in more detail. Half the usual number of gradient steps with twice the 
usuaFspacing in k-space were applied, and the RF response was recorded in each 
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of three coils. Given the gradient phase encoding regimen just described, each of 
these RF responses corresponded to an image of half the desired J field of view. 
The coil signals were then corhbined into two 'sets with different weightings to 
produce synthetic composite signals corresponding * to" tKe 'zero* 11 and first 
5 harmonics;' and the tivo synthesized fespbr^es : were interleaved to fill alternating 
lines' of the" k-space matrix/^ Founer^ tfansforrned in a 

conventional 1 way t6 form the image. The resulting image had Tull resolution 
over the full field of view. As shown in FIGURE 5, line (A) depicts* the 
acquisition of signal in each compqnenLcoiL A>representatiye signal point P is 

10 picked up in each coil, and thus contributes to the line in k-space formed with 
each. coil signal. Line (B) shows the processing performed in the prototype 
.described above, according to which the. coil signals, are combined into a first set 
S 0 and a second set S r Avhich each: form different lines (even and odd as shown, 
corresponding to one gradient step separation) of, the ; k-space matrix. These are 

15 combinedat line (C) into the full, matrix. : ; f c ; 

FIGURE ? 5A illustrates corresponding images and signals, Line (A) of 
that FIGURE sho>\vs the half field-of-view 'images a, b, c reconstructed from 
separate coil signals, of which the signals themselves are shown in iine (B). Line 
tQ shows : tfte two spatidl harmonic "signals and iine (D) the full interleaved 
20 signal, or k-space matrix. Line (E) illustrates the image reconstructed from the 
signals of line (C)V The middle three Stages (B), (C) and M (D) correspond directly 
to the steps shown schematically iil FIGURE '5" ^ - • ' > • f - 

A ^ ■ FIGURES 6A, 6B and -7 A; 7B, respectively^ show phantoiri and in vivo 
images formed by the three-coil prototype embodiment using the SMASH 

25 imaging technique. FIGURE 6A shows? the; proton-density image of a water- 
filled phantom, obtained using a TSE sequence (described further below) with a 
total conventional acquisition time of 10 seconds." FIGURE 6B shows the 
equivalent image obtained in 5 seconds with the SMASH processing of the 
present invention. Apart from certain residual artifacts, the images have an 

30 identical appearance. 

Similar results were obtained in. vivo. FIGURES 7 A and. 7B shqw a 
reference and SMASH-reconstructed coronal image through the brain of a 
healthy .adult volunteer, the images being acquired in; 71 ..and ., 35 seconds, 
-respectively. . All these images were acquired using commercial hardware and a 
35 convenient, pulse sequence with high spatial resolution \and -good SNR; a fifty 
percent acquisition time reduction was achieved. The acquisition time savings 
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described here have also been shown, to apply, to other commercial machines and 
pulse sequences, including the. fastest MR , imaging machines and sequences. 
Some residual foldover, artifacts are.present in the. SMASH. reconstructed image, 
due to, imperfections in the . cprnposite spatial .harmonic sensitivity profiles. 
5 However, the .prototype, used.an existing .poil array designed. for other purposes, 
and improvements., both in coil ;l desigru and, in the, accuracy of coil sensitivity 
mapping and j spatial harmonic generation, are .expected to minimize these 
artifacts.- , * s > r > ; - . v . j(/: - : - \ ; . ./ * tf * 

'•• The'raw data* for -the foregoing images were generated dn a commercial 
10 1.5 Tesla whole 'body clinical MR imager, which was a Philips NT, made by 
Philips Medical * Systems of Best, Netherlands. For the phantom images, a 
standard circular phantom of 200 miti diameter and varied internal structure was 
used.^ A» prototype "coil array on hand in applicant's laboratory was used as the 
SMASH coil array. The 1 equivalent circuit for this coil consisted of three 
1 5 rectangular component coils arranged in a linear array with a slight overlap in the 
right-left direction, similar to the array shown schematically in FIGURE 2B. As 
these coils were designed for another purpose, details of the coil circuitry, not 
relevant : here resulted in more complicated sensitivity profiles than would Jbe 
expected from the simple geometry . described above. It was, these sensitivity 
20 variations which were primarily , responsible- for the residual, foldoyer artifacts 
noted in FIGURES 6B and 7B.,-A simpler, coil designed : with spatial- harmonic 
generation in mind will circumvent-these artifacts. . . r . 

; Further, details of the prototype imaging proces? and equ r ipmerit are given 
.below., ^ , _ v,-,. ,-, .1, v <; ' 

25 v * * Image planning and data 1 collection ^ 1 *'V vi \;> :.».- 

The phantom images of FIGURE 6A, B were generated as fol?6ws. First 
with the phantom centered over the coil array, data for the reference image were 
acquired in a 6 mm thick coronal slice parallel to and 80 mm above the plane of 
the array, using a turbo spin echo pulse sequence with five echoes per excitation. 

30 The field-of-view (FOV) was 200 mm, centered on the phantom, and matrix size 
was 256 ; x 256. Phase encoding was performed in the^right-left direction (i.e. in 
the direction of the coil array); with a single signal average. Data from each of 
the three component coil channels was acquired simultaneously and,' in the 
prototype, were stored separately for later processing. Acquisition time was 

35 measured at 10~seconds. ' Next, a second coronal slice using the same technique 
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and imaging parameters was taken to serve as a measure of component coil 
sensitivity! It was acquired 12 mm above the first one, in a region of uniform 
spin density in the phantom; r Then, a third image at the kame level as the first 
was obtained in half the time using twice the; phase-encode step, and hfence half 
5 the field of view, in the right-left direction. Matrix : size was now 256 x 128. 
'* Acquisition time Was 5 seconds, exactly half thb iime taken for the 1 first image. 

For the in vivo images of FIGURE 7 A, B, the volunteer was positioned 
with his head above the coil array, and images were taken in the same plane and 
with the same parameters as for the phantom images, ekc6pt that eight signal 
10 averages with a slice thickness of 10 mm were used to improve SNR. The coil v 
sensitivity image from the phantom was also used as a sensitivity reference for 
the in vivo images. : u/ 

Determination! of coil sensitivity and optimal component coiO weights 

As described above in connection wjth FIGURE 4, coil weightings were * 
15 determined by fitting actual coil sensitivity data to, the desired spatial harmonic 
sensitivities. To do this, the sensitivity function along the right-left diameter of 
the phantom was extracted from the component coil sensitivity reference images. 
The phantom spans the entire extent of each image, and since the sensitivity 
' ' reference image plane intersects the phantom in a region of uniform spin density, 
20 these intensity profiles correspond precisely to the complex sensitivity functions 1 
' * 6f the coils along the diameter. 1 The' sensitivity of each component surface coil 
falls off rnohotonically with distance from the coil, and' tht phase shift grows 
with distance from the coil center; Using a gradient-descent fitting routine, these 
' sensitivity functions were 1 iterativdy fitted to twb targe* spatial harmonics: to the 
25 ! zero^frequency harmonic corresponding to uniforiii Verisitivity/ and to the first 
- harmonic having a modulation wavelength equaUo theTield-of-vieW of 200 mm. 
The fit, shown in panels C and 'D of FIGURE 4A, demonstrates that veiy close 
fits to the target harmonics may be achieved by component coil weighting, even 
: with only three coils 1 . A slight residual ripple was present in the zero th harmonic 
30 ' fit and this is primarily responsible for imperfections iri the SMASH image - 
reconstructions. ~ ' 

Image reconstruction 

The reference images in FIGURES 6A AND 7 A - were formed by 
* combining, the component coil reference images pixel-by-pixel as the square foot 
35 of the sum of square magnitudes. *• _~ . , > , . 1 ■ \ • 
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FIGURE 5 A shows intermediate stages in SMASH reconstruction of the 
phantom image. The, procedure for the in. vivo image was, identical. Using 
weights from the - iterative fit, the three .component, coil signals (B) representing 
half-time, .halfrFOV aliased images (A) - : were combined into two composite 
5 :. signal sets, one, for the zero** 1 spatial; harmonic and , one , for the first spatial % 
harmonic (C).» .Finally,, the two. composite signal .data sets were interleaved to 
form a data matrix of size 256 x 256 (D), This matrix was Fourier transformed to 
yield the reconstructed image (E), which is also shown in FIGURE 6B. 

In general, the foregoing SMASH imaging technique may be applied with 
10 a great many.of the, known pulse sequences or spin conditioning; techniques, and 
will in general share the advantages, of the underlying sequential imaging 
methods which are used to collect partial £-space information. , As long as 
suitable spatial harmonics may be generated with a coil array, the additional 
acquisition time savings afforded by SMASH reconstruction involves no 
15 significant sacrifice in -"resolution or SNR. " This -contrasts ' markedly with the 
trade-off in SNR or resolution that characterizes many existing approaches to fast 
imaging, such as low flip-angle sequences:; ' » ■ . . . . = 

Signal-to-noise ratio 

Visual inspection of the, phantom, images in FIGURES 6A, 6B reveals a 
20 slight degradation in SNR in the SMASH reconstructed image as r compared with \ 
.the reference image. Part of this apparent loss may be traced to. the., residual 
. foldover artifacts, in which some of. the intensity of the primary, image, is "stolen" 
. by aliased ghosts. There : is .another noteworthy difference between the ; SNR 
xf profiles of the two images, however. The reference image . in.,FIGURE t 6A .was 
25 , _ generated using a, sum-of-squares combination of component coiljmages, ^s is ; r ; 
. described. jn PB. Roemer, W.A. Edelstein, C.E. Hayes, S.P. Souza, and O.M. 
Mueller, The NMR phased array . Magn. Resort. .Med, 16, 192-225 (1990) and as 
has become a standard. practice in phased array imaging to improve SNR. jThis 
^combination algorithm yields an essentially constant noise profile t across .the 
30 . image, with.varying signal (and hence varying SNR) that is enhanced in regions 
of significant overlap between component coil sensitivities. SMASH 
reconstruction, on the other hand, results in a consXani-signal image with a 
varying noise profile. The linear combinations of component coil signals in 
SMASH are explicitly designed to produce a homogeneous composite signal. 
35; This is most obvious in the case of the zero 1 * 1 harmonic combination, which 
produces a flat net sensitivity profile C tot across the image plane/ While 4hc - 
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higher- harmonics do involve significant spatial sensitivity variations, their 
profiles are all complex- exponentials of •unit -modulus; and they do not lead to 
any intensity variations in the absolute magnitude image: (In other words, any 
apparent loss of sensitivity ' in the real channel is precisely compensated by a 
5 'sensitivity gain an the imaginary -channel.) SMASH reconstructions, therefore, 
' do awaj"with ihtensityrpeaks. in- regions of component' coil overlap in favor of a 
: spatially * homogeneous Mm'ajge profile; Thus, sensitivity-dependent linear 
combinations may be * used independently ^ of - SMASH ? *as ^a method of 
homogeneity correction^ A .homogenejty-cflrreqted version of the image in 
10 FIGURE 6 A may be generated by, linear combination, of .component coil 
reference images, using the weights calculated for the zero**? harmonic profile in 
FIGURE 4A. Roemer et. al. have described an alternative combination 
. algorithm for producing constant-signal images from phased arrays, and SMASH 
images produced with accurate spatial harmonics will have ; similar intensity 
15 .profiles. .. . .... . , 

• Such considerations of -image intensity profiles arfe addressed primarily at 
the numerator of the signal-to-noise ratio. The behavior of noise in SMASH 
images is complicated by the presence of noise correlations among the vkriously 
combined data sets that enter into a SMASH reconstruction. Noise from each of 
20 the individual, coils is replicatqd (a}beit :: >yith . different weights) in each of the 
.different harmonic composite signals, ,and,this.niight raise a concern that some 
portion of tlie noise voltages will .add^coherently. If the , addition were purely 
coherent, i.e. if nqise were to enter into each, of, the composite signals with 
identical .vyeights, then the ,ims v ? poise^ value \y014d exceed the value for 
25 ^conventional full-time acquisitions t by > the., square- root of the number of 
composites. .Jtn other words, the upper limit pf SNR loss due, to noise correlations 
,is the square root of the time sayings factor. , In practice, however, the signal 
combinations arc formed, with distinct sets of complex weights, which -reduces 
, noise? correlations substantially, below this^theoreticaflimit; t , 

30 Dependence on component coil sensitivkies • ' fc • 

In the above described embodiment, the SMASH imaging technique 
depends upon the measurement arid manipulation of "coil' sensitivities for its 
spatial harmonic generation. Errors in spatial harmonics may lead to aliasing 
- artifacts and Sl\fR degradation (for example due to mismatch arid mixihg between 
35 the 'spatial frequency components of the image): This raises severar practical 
'issues. First, we-must be able to estimate coil sensitivities with some reasonable 
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degree of accuracy. .Sensitivity maps from phantom images were used for the 
, SMASH reconstructions presented here. As noted above,- numerical simulations 
of coil sensitivity may, also be -used, . if these simulations, have been reliably 
calibrated. The invention has also implemented use of other techniques for direct 
5 sensitiyity estimation, .including* in vivo, estimations using .fast low .resolution 
scans at the time of iw, vivo imaging. Lqwrpass-filtering of images could also be 
used to provide sensitiyity maps. These techniques may ? be .^automated for 
calibrating, normalizing and deriving weights /overeating thecomposite.signals. 

Another detail to be addressed * is that Bq and - B ) magnetic field 
10 inhomogeneities m&y distort the true coil sensitivity profiles and interfere with 
spatial harmonfc generation. Spin echo sequences, like the TSE sequence used in 
this work, have the advantage of refocusing static field inhomogeneities. 
'Gradient echo sequences have been' found to perform just as well. The weight- 
finding algorithm used by applicant automatically compensated for some degree 
15 of residual field inhomogeneity and generated accurate spatial harmonics even in 
the presence of systematic phase or intensity errors.. In fact the technique has 
been found to-be directly applicable not only to spin echo sequences^such as TSE 
but also to several variants of gradient -echo, sequences including FISP and 
\ standard GRE. ■ . » \ . 1 ir r - .«■>.- 

20 Finally; the geometry of the coil array will place certain limitations on the 

field of view (FOV), the positibn, and the angulation of image planes suitable for 
SMASH reconstruction/ Surface coil sensitivity profiles vary with distanbe from 
the Coilc'enter, and while simulations indicate that a wide range of image-plane 
geometries are compatible with-SM ASH reconstruction; the reconstiruction mfay 

25 begin to fail at large distances" and angles, where sensitivity functions become 
broad and asymmetric. These constraints-will be relaxed iri the presence of larger 
numbers of component coils whose sensitivities provide good coverage of the 
imaging volume. The design of RF coilarrays specifically tailored fof'accurate 
and flexible spatial harmonic generation' will* help to overcome the limitations 

30 outlined above. - . r f 

In particular, the invention contemplates the use of coil arrays with 
multiple component coils extending in more than . one linear direction. A two- 
dimensional array, such, as .an, N xM rectangular array, will allow generation of 
spatial harmonics along multiple directions and will relax constraints on image 
35 plane position and angulation. Wraparound arrays are also contemplated to allow 
spatial harmonic: generation in a f plane transverse tothe body plane. A coil 
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arrangement with a surface array on the top of thew body and : another surface 
array on the botto'm, possibly withsome linear offset with respect to one another, 
-will allow fine tuning of spatial harmonics in a plane between the two and will 
' have the added advantage of iricreasiing overall SNR in- siich- a' plane/ An 
5 extended grid coil grid coil may also be designed'withf a sensitivity profile more 
clpsely^ ( apprqximating a . sinusoid, AH --such; multiple-coil configurations will 
. allow improvement in the generation of spatial hannonics, as.vvell as, allo\ying a 
, greater; number and- variety of harmonics ,to : be . generated. As in ; the prototype 

- embodiment described, , the component xoi Is in 'these arrays may be overlapped 
10 and may be output to low input impedance preamplifiers as necessary,- to 

- minimize inductive coupling. It is -also possible that the simple, expedient of 
transforming, the coil, sensitivities to achieve the jspatial harmonic reconstruction 
procedure of applicant's SMAg.H invention mayvbe applied, in large arrays of 
coilsof the % sort proposed by MR; Hutchinson and Raff in 1988 and: described 

1 5 . further byKwiat, Einav and Navqn, supra..: : - 

: It should be observed that, although recent studies have shown that the 
sensitivity advantages of phased arrays peak at coif numbers less than 10; the 
advantages of faster imaging with SMASH imaging may be -further enhanced by 
'arrays wittf : larger : numbers of coils arid with new" coil geometries/ As one may 
20 intuit from FIGURES 1 3 ahd T 4A, the higher the ffbqiiericy ; of the spatial harmonic 
desiPed/ thd larger the ' number of component 1 coils that will be required for its 
1 generation; This is'in part : because, at least for simple Component coil Shapes, the 
Sharpest features bf the" net sensitivity fuhctioh " carl orily be on tfie order of the 
cbmpbnent cbif dimension. As rioted' above, the' extrk sensitivity information 
25 frorii additibrial coils can also be'used fortHe-fine tuning of spatial harrrioriids ^id 
5 the improverifent'of image quality. J r * ' : - yrt * v - • ^ 

. ? The SMASH technique as applied to a conventional NMR apparatus 
described above partially replaces gradient phase encoding by a spatial encoding 
procedure tied to the detection coils, in which some of the spatiai modulations 

30 that distinguish different phase-encoding lines are generated by combining 
signals from multiple coils arrayed above or around the imaging volume. This 
shift of responsibility from gradient geometry to coil geometry, and from the spin 
preparation stage to the stage of signal detection arid combination, allows for 
simultaneous acquisition of multiple lines of £-space. Enhancement 6f image 

35 . speed, resolution or ; field of view by a factor of two,, fiver ten or more may be 
expected with appropriate coil arrays generating sufficient harmonics. 
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Apart from its promise for progressively faster acquisitions, SMASH has 
a, number of practical t advantages as a fast imaging- scheme. . As a partially 
parallel acquisition strategy, it can be, combined with most existing sequential 
fast imaging techniques; for multiplicative time savings^ No. special hardware is 
5 required, other than an appropriate coil array ... A > t . .... 

One embodiment* of thie invention includes a coil 'array together with a 
digital- signal processor- configured to , combine the outputs of the component 
coils with appropriate weights as' described above, and produce -two or more 
output signals; each of which represents a composite spatial harmonic as 

1 0 described above. The array with processor may then be directly substituted for a 
conventional receiving coil, with the difference that it produces two or more lines 
of data for eabhfspin preparation cycle of the MRI apparatus. As such, the plug- 
in coil unit ' may operate with older MRI devices to directly enhance the 
acquisition time, field of view; or resolution by a factor of two or more, without 

15 changing the expensive magnet and other spin preparation hardware of the 
device. Commercially available or . home-built >MR phased arrays may well 
^suffice to yield significant time savings, and SMASH may be performed on 
. machines not. equipped with much more costly magnetic enhancements such as 
. EPI gradient systems. Combination of the component . coil signals , may. ; be 

20 .performed after, the fact, allqwing for. a -wide, range, of postprocessing . steps, 
. .including, fine tuning. of spatial harmonics, adaptive artifact, correction,.- or even 
. nonrFourier (e.g. wavelet) encoding and reconstruction^ Furthermore. t iL should 
be . noted that the inventipn further contemplates, that both for constructing .the 
basic component ..coil transforms .or weights and for any. timinjg or adaptive 

25 adjustments, the. invention, may -include ; a dighal -signal .processor ; -pr. neural 
network processor in addition to the normal . complement of, control , and 
processing assemblies to evaluate the responses or coil signals and carry out the 
signal combination into spatial harmonics in real time. Even in the simplest 
embodiments, spatial harmonic reconstructions automatically yield a 

30 homogeneous intensity profile for" the reconstructed image, and this may be 
advantageous for some imaging applications. 

The time gained in a SMASH acquisition may be used to collect extra 
data for better spatial resolution and/or SNR, or else faster acquisitions may be 
used to eliminate motion x artifactsirom mobile structures in the field of view . , 

35 The invention \ has been described above with reference to : a prototype 

embodiment which . produces, signals corresponding to spatial harmonics that 
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yield additional k-space lines equivalent to the conventional gradient steps in a 
Fourier transform imaging system and a simple coil array and processing 
assembly for enhancing image speed, quality or field. However, the method is 
broadly applicable to coils and processors for "producing various spatial 
5 representations, including wavelets or other non-Fourier spatial representations to 

' 1 'replace, supplement' or enhance- MR imaging processes 3 ahd systems. The 
invention being thus described,TurtHef variations arid m occur to 

those skilled in the art, and all such variations and modifications are considered 
to be wiWiiri the scope ; of the invention and its equivalents, as set forth in the 

10 claims appended hereto. ! : * * » 

, ■* j ... v • .... .» 

What is claimed is: 
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Claims * ........ 

; 1 . f A method pf imaging a region of a body by MR by capturing- MRI data during 
multiple cycles, of NMR spin conditioning and measurement, said method comprising 
5 the steps. of: «. • . . .... . f .. : - ; - - •>...,-. ' . . - ^ 

i)^ ; , during one spin conditioning cycle receiving an.NlylR RF. measurement 
signal response simultaneously, in each of plural surface coils haying different spatial 
sensitivities,. ... v ,.,~ - -. ; : .... 

, ii) ... .combining said NMR RE measurement, signal responses received in the 

1 0 plural coils into plural sets of composite NMR signals, each set of composite NMR 
signals containing different spatial information, 

iii) repeating said steps of receiving and combining during additional spin 
conditioning cycles to form a data set augmented by at least one set of said composite 
NMR signals, and 
1 5 iv) forming an MR image from said data set. 

2. A method according to claim 1 , wherein said steps of receiving and combining 
are carried out to form images of a given field of view and spatial resolution while 
reducing at least one of number or duration of spin conditioning cycles and thereby 

20 reducing image acquisition time. 

3. A method according to claim 1, wherein said steps of receiving and combining 
are carried out to form images in which at least one of, field of view, spatial resolution, 
and acquisition speed is enhanced. 

25 

4. A method according to claim 1, wherein the step of combining combines the 
measurement signal responses in linear combinations. 

5. A method according to claim 4, wherein the linear combinations are performed 
30 using complex coefficients. 



6. A method according to claim 1, claim 4, or claim 5, wherein the step of 
combining signal responses yields composite spatial sensitivities which approximate 
spatial harmonics of an imaged field of view. 

35 
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7. A method according to claim 6, in which the spatial harmonic composite signals 
are substituted for data lines corresponding to omitted increments in magnetic field 
gradients. .. . . * 

'5 8. A method according to' claim 6; wherein the step of receiving in surface coils' 

includes receiving in at least one of a two dimensional coil array* an Nx Mcoil array, a 
wrap around coil array and an extended grid coil array. 

9. A method according to claim 1 , further comprising the step of sensitivity 

1 0 mapping said plural coils and determining a transform for combining the received NMR 
RF measurement signal responses into said sets of composite NMR signals. 

10. A method according to claim 6, wherein the cycles of spin conditioning and 
measurement are selected from among spin echo, gradient echo and BURST style 

15 imaging techniques. - . . - 

11. A method according to claim 1 , wherein the step of combining into sets of 
composite NMR signal responses combines into sets of signals corresponding to non- 
Fourier spatial distributions, and the step of processing applies ai corresponding non- 
20 Fourier transformation to form the MR image. ' ,; ' " 1 J 

12. A method according to claim 1 ^ wherein the step of combining into sets of 
composite NMR signal responses combines into wavelets forming sets' of spatially 
Orthogonal data, knd the step of processing applies a wavelet transform to form the 

25 image. * ' '' ; 

' v l3. A magnetic resonance imaging XMRI) apparatus'of the type having a magnet for 
establishing a magnetic field in a region means for conditioning nucleiar spins in the 
region and means for collecting an RF signal from the conditioned spins, said magnet 
30 and said means being operated in repetitive cycles to collect spin data iand construct ah 

image of material in said region; wherein said means for collecting an RF signal includes 

a plurality of surface coils forming an array, and each having a different 
spatial sensitivity, arid ' * ~ * ' ' 

means for combining signals from plural individual coils of said array 
35 with different weights into a plurality of sets of signals corresponding to orthogonal 
spatial representations of spin data from said region. 
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14. A magnetic resonance imaging (MRI) apparatus according to claim 13, wherein 
said means for combining combines said signals in flow-through fashion as said spin 
data are collected.. 

5 15. H . A.magnetic resonance imaging (MRI) apparatus according,to claim 13, wherein 
said spatial representations are spatial harmonics.. . , f .... .... . . , 

16. A magnetic resonance imaging (MRI) apparatus according to claim 13, wherein 
said spatial representations include a non-Fjourier representation.,,: . ^ 

10 . .... 

17. A magnetic resonance imaging (MRI) apparatus according to claim 16, wherein 
said non-Fourier representations are wavelets. 

18. A magnetic resonance imaging (MRI) apparatus according to claim H, further 

1 5 comprising means for determining said weights, and wherein said means for combining 
combines in accordance with weights determined by said means for determining. 

19. . A magnetic resonance imaging (MRI). apparatus according. to claim 18, wherein 
said means for determining includes a stored table of weights for each of $ plurality of 

20 imaging regions for forming said spatial representations. , .. - ;r , . . , r 

. . 20. , A magnetic resonance, iniaging (MRI).apparatus according to claim 1 3 , wherein 
said plurality of surface coils form an array having ^extent of approximately twenty to 
about fifty centimeters, and. individual coils of the array have an extent each have an 
25 extent under about twenty centimeters. 4 . r : 

21 . A magnetic resonance imaging (MRI) apparatus according to claim 13, wherein 
_ said plurality of surface, coils form a linear array. 

30 22. ; In a magnetic resonance imaging (MRI) apparatus of the type having a magnet 
for establishing a magnetic field in a region, means for conditioning nuclear spins in the 
region and means for collecting an RF signal from the conditioned spins, said magnet 
and said means being operated in repetitive cycles to collect spin data and construct an 
image of material in. said region, the improvement wherein said means for collecting an 

35 RF signal includes , .... 

a plurality of surface coils forming an array, and each haying a different 
spatial sensitivity, and 
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means for combining signals from plural individual coils of said array 
with different weights into a plurality of sets of signals corresponding to orthogonal 
spatial representations of spin data from said region. ; 



5 ; * 
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